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ABSTRACT 


The  Green's  function  approach  to  the  temperature -dependent  many 
body  problem  is  reviewed  in  some  detail  (Chapter  II).  A  compact 
formulation  of  the  problem  is  proposed  which  emphasizes  the  well-known 
parallelism  between  the  roles  of  the  time  and  a  certain  temperature  - 
variable  in  the  temperat  ure -dependent  many  body  problem.  This 
formulation  is  then  used  to  derive  in  a  new  way  some  formulae 
previously  established  in  the  literature  (Chapter  III).  The  energy- 
momentum  representations  for  the  Green's  functions  are  then  obtained, 
and  the  parallelism  which  was  mentioned  above  is  discussed  A  result 
of  this  discussion  is  to  show  that  the  parallelism,  though  inexact,  is  far- 
reaching.  In  this  connection  it  is  shown  that  the  Green's  function 
pertaining  to  the  temperature  variable  can  be  split  into  two  parts  which 
can  be  associated  unambiguously  w  th  the  retarded  and  advanced  parts 
of  Green  s  functions  pertaining  to  the  time  (Chapter  IV).  In  the  last 
chapter  the  main  results  of  this  thesis  are  discussed  in  the  context  of 
the  recent  literature  and  some  suggestions  are  made  for  further 
investigations  . 
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CHAPTER  I:  INTRODUCTION 

In  dealing  with  the  many  body  problem  (MBP)  one  is  concerned 
with  an  ideal  zed  physical  situation,  viz.,  one  investigates  the 
properties  of  a  system  composed  of  N  particles  in  a  volume  V  where 
both  N  and  V  are  infinite  but  their  ratio  (N/V)  is  finite.  It  is  clear 
from  the  outset  that  for  a  system  with  a  large  number  of  interacting 
particles  one  cannot  hope  to  solve  the  Schrodinger  equation 

H f  = 

exactly,  i.  e.  ,  one  cannot  expect  to  get  a  solution  for  \j/  in  closed 
form.  Furthermore,  we  are  not  really  interested  in  the  detailed 
solution  ;  what  is  of  interest  to  us  in  the  many  body  problem  is  to 
explain  the  gross  properties  of  the  system  under  consideration  in  terms 
of  some  fundamental  properties  of  the  particles  which  make  up  the 
system.  To  be  specific  the  fundamental  properties  that  are  of  interest 
here  are  the  mutual  interactions  among  particles  and  the  statistics 
obeyed  by  the  particles. 

As  a  result  of  the  mathematical  complexity  of  the  problem  it  is 
customary  to  use  models  to  describe  a  system  composed  of  a  large 
number  of  interacting  particles.  Thus,  the  fission  of  nuclei  is  pictured 
as  the  fission  of  a  drop  of  liquid.  Another  process  for  which  a  liquid 
drop  model  of  nuclei  is  useful  is  the  photoemission  of  nucleons  from 
heavy  nuclei.  As  is  to  be  expected,  this  model  has  its  limitations;  for 
example,  the  photoemission  of  high  energy  nucleons  from  heavy  nuclei 
cannot  be  explained  in  terms  of  the  liquid  drop  model. 


The  purpose  of  a  satisfactory  many  body  theory  in  the  case  of  zero 
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temperature  is  to  obtain  a  convenient  expression  for  the  ground  stale 

of  the  system  in  terms  of  the  fundamental  properties  of  the  particles 

which  were  specified  above.  Actually,  when  one  wishes  to  calculate 

wi  h  this  expression  the  characteristics  of  a  given  process,  one  must 

resort  to  approximations.  These  approximations  may  be  different  for 

different  processes.  In  executing  these  approximations,  one  should  be 

in  a  position  to  determine  the  region  of  validity  of  a  given  approximation. 

Since  each  approximation  corresponds  to  a  model,  one  can  estimate  from 

the  outset  the  validity  of  a  model  for  a  particular  process.  For  example, 

the  Brueckner  theory  was  shown  by  Goldstone  (19  57)  to  take  into  account 

a  certain  class  of  terms  in  the  general  expression  for  the  ground  state 

energy  of  nuclear  matter.  These  terms  were  shown  by  Brueckner  to 

justify  the  single  particle  model  for  such  a  system.  Here  we  have  an 

example  of  a  many  body  theory  which  justifies  a  known  model  for  nuclear 

matter;  moreover,  in  this  example  the  theory  actually  enriches  the 

potential 

model  inasmuch  as  it  shows  that  the  appropriate  single  particleAis  velocity 
de  pendent . 

The  considerable  progress  in  many  body  theory  which  has  occurred 
in  recent  years  is  due  to  a  large  exvent  to  the  realzation  that  the  arsenal 
of  field  theory  can  be  used  in  an  attack  on  the  MBP.  This  realization 
is  usually  attributed  to  Goldstone  (1957)  who  used  the  techniques  of  field 
theory  employing  the  so-called  interaction  picture.  However,  when 
applied  to  the  MBP,  field  theory  techniques  employing  the  Green's 
function  formalism  possess  a  distinct  advantage.  This  is  primarily  for 
two  reasons:  fa)  In  the  Green's  function  formalism  one  considers  only 


the  so-called  linked  cluster  diagrams.  These  diagrams  correspond  to 
processes  wherein  particles  are  excited  from  the  ground  s*ate,  and 
after  going  through  a  variety  of  intermediate  sta'es  excluding  he  ground 
s'ate,  return  to  the  ground  state.  In  the  alternative  apprpach  using  the 
interaction  picture,  one  obtains  terms  which  correspond  to  processes 
wherein  particles  return  to  the  ground  s'ate  before  the  end  of  the  process 
("unlinked  clusters").  Although  we  expect  that  the  terms  corresponding 
to  unlinked  clusters  can  be  eliminated,  the  procedure  of  proving  this 
elimination  is  not  rigorous  (Goldstone,  1957;  Bloch,  1958)  (b)  The 

Green  s  function  formulation  is  concise  and  elegant .  For  example,  in 
this  dissertation  the  whole  problem  has  been  formulated  without  recourse 
to  diagramatic  techniques.  At  this  juncture  we  would  like  to  uote  the 
motivation  given  by  Klein  and  Prange  (1958)  for  introducing  the  Green's 
function  formalism  to  the  MBP;  "The  formal  development  of  the 
subject  (MBP)  has  been  marked  by  an  increasing  reliance  on  the  method 
of  second  quantization,  not  at  all  surprising  in  view  of  the  fact  that  we  are 
dealing  from  one  standpoint  with  a  class  of  particularly  simple  field 
theories.  What  has  seemed  surprising  to  the  authors  is  that  detailed 
treatments  so  far  published  have  for  the  most  part  been  either  special 
to  the  problem  at  hand  or  else  have  taken  over  general  techniques  of 
quantum  field  theory  which  date  from  the  'early  part'  of  the  last  decade 
of  profound  development  of  the  methodology.  " 

In  this  thesis  we  concern  ourselves  with  the  temperature -dependent 
many  body  problem  (T  -  dependent  MBP).  In  the  T -dependent  MBP  we 
are  interested  in  the  function  which  gives  the  density  of 
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levels  as  a  function  of  the  energy  and  the  number  of  particles  per  unit 
volume.  A  related  function,  with  which  we  will  deal  directly,  is  the 
grand  partition  function  Z.  This  latter  function  is  defined  by 


where 


Z  =  77.  , 

Z.r  =  sxpbH +*N]. 


Here  H  is  the  Hamiltonian  of  the  system,  N  is  the  particle  number 
operator,  jS  -  k  /  (T  being  the  temperature  of  the  system  and  k 
the  familiar  Boltzmann  constant)  and  o(  is  related  to  the  chemical 
potential  yH  by  yM  ~  °^/ fi>  ■  The  grand  partition  function  may  be 
considered  as  the  fundamental  function  of  the  T  -dependent  MBP  since, 
as  is  well  known,  all  the  thermodynamic  properties  of  a  many  body 
system  can  be  derived  from  Z. 


The  operator  /iop  satisfies  Bloch's  equation: 

*LKcr  -  -  N  ?  . 

This  equation  is  very  similar  inform  to  the  Schrodinger  equation.  This 
similarity  has  been  utilized  by  various  authors  to  develop  an  expression' 
for  Z  in  terms  of  the  mutual  interaction  potentials  of  the  particles, 
e  g.,  Bloch  and  DeDomiricis  (1958)  and  Luttinger  and  Ward  (I960) 

(these  papers  will  be  referred  to  as  BDD  and  LW  respectively).  This 


procedure  is  mathematically  similar  to  the  method  in  field  theory  which 
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is  known  as  the  interaction  picture  and  which  was  mentioned  above  in 
connection  with  the  zero  temperature  case.  In  this  approach  Bloch  and 
DeDominicis  used  an  exact  generalization  of  Wick  contractions  (Wick,  1950) 
applicable  to  the  T-dependent  system  with  the  many  body  system 
considered  in  a  grand  ensemble.  BDD’s  expression  for  the  energy  of 
the  system  bears  strong  resemblance  to  Goldstone  's  expression 
(Goldst  one,  1957)  for  the  ground  state  energy  of  the  many  fermion  system. 
(In  this  dissertation  we  deal  exclusively  with  fermions.)  However,  it 
was  pointed  out  by  Kohn  and  Luttinger  (I960)  that  in  actual  fact  the  BDD 
expression  does  not  reduce  to  that  of  Goldstone  in  the  T  -*■  O  limit  , 

Kohn  and  Luttinger  showed  that  in  taking  the  T-**0  limit  in  the  BDD 
expression,  one  gets  ,  in  addition  to  the  Goldstone  terms,  the  so-called 
anomalous  terms. 

Luttinger  and  Ward  (I960)  and  Balian,  Bloch  and  DeDominicis  (1961) 
showed  that  for  a  central  two  body  interaction  (i.e.  ]/  (  }  X-  )  = 

\/  (  l  X- ~  Xj  |J  where  V  (  % j)  is  the  interaction  potential 

between  the  particles  ,  i  and  j  )  the  two  expressions  -  Goldstone 's  and 
BDD' s  -  can  be  brought  into  coincidence  ,  i.e.  the  anomalous  terms 
can  be  removed  if  one  allows  for  the  possibility  that  the  chemical 
potential  is  different  in  the  two  cases. 

In  the  LW  presentation  extensive  use  is  made  of  the  "single  particle 
propagator.  "  In  their  context  this  propagator  traces  what  one  might 
call  "the  propagation  of  the  particle  in  temperature".  In  this  formalism 


. 
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a  variable  which  is  related  to  the  temperature  plays,  formally,  the  role  of  time 
and  the  corresponding  propagator  is  closely  related  to  BDD  contractions. 

Matzubara  (19  5C)  was  the  first  to  show  that  it  is  possible  to  define 
a  temperature  dependent  Green's  function  which  possesses  formally  the 
same  mathematical  properties  as  the  Green  s  function  that  wa6  introduced 
in  field  theory  by  Schwinger  (1951).  A  stimulating  generalization  of  the 
Matz  Tbara  approach  was  worked  out  by  Fradkin  (1959)  who  considered 
grand  ensemble  averages.  Fradkin  showed  that,in  his  formulation,  the 
time -de  pendent  Green's  function  for  a  temperature-dependent  problem  is 
most  naturally  defined  as  the  analytic  continuation  of  an  appropriate 
T-dependent  Green  s  function 

Apart  from  their  aesthetic  appeal,  and  the  possibility  that  future 
developments  may  enhance  their  value,  one  might  question  whether  the 
method  of  the  temperature  Green's  functionsis  useful  in  its  own  right. 

That  it  is  useful  can  be  seen  from  the  following  argument.  The 
expression  for  the  single  particle  propagator  in  time  (referred  to  as  the 
"time  Green's  function")  is  known  to  have  a  direct  physical  meaning. 

This  can  be  stated  most  succinctly  in  the  following  way:  the  poles  of  the 

time  Green  s  function  in  the  energy -momentum  representation  g  ve  he 
energy  spectrum  of  the  system.  However,  evaluation  of  the  expressions 
obta  ned,  when  one  attempts  to  solve  the  propaga  or -in-time  equation 
directly,  appears  difficult  if  not  impossible  (Landau,  1958).  Thus  it  is 
very  useful  to  obtain  via  some  other  procedure  expressions  which  seem 


' 
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simpler  and  yet  are  capable  of  yielding  the  same  information.  For  this 
reason  the  development  in  the  temperature  variable  is  of  interest.  It 
was  shown  by  Abrikosov,  Gorkov  and  Dzyaloshinski j  (1959)  (this  paper 
will  be  referred  to  as  AGD)  that  the  retarded  time  Green's  function 
or  propagator  is  simply  related  to  the  temperature  Green's  function  in 
the  energy -momentum  representation. 

The  present  investigation  was  motivated  by  the  striking  similari  ty 
between  the  roles  of  the  temperature  and  time  variables  in  the  many  body 
problem.  A  central  purpose  of  this  work  has  been  to  investigate  further 
th  s  interesting  similarity.  In  this  connection  we  have  investiga  ed,  in 
energy -momentum  representation,  "retarded  and  advanced"  temperature 
Green's  functions  in  analogy  with  the  time  Green's  function  along  the 
lines  given  in  Fradkin’s  paper  (Fradkin,  1959)  To  this  end  a  new  formal 
scheme  was  devised  whereby  the  parallelism  of  time  and  temperature  is 
manifested  in  a  clear  way.  For  example,  the  formalism  contains  an 
"equation  of  motion"  for  the  Green's  function  in  the  temperature  variable. 
The  analogy  is  further  enhanced  in  Chapter  III  where  it  is  shown  amongst 
other  things,  that  the  temperature  variable  can  formally  be  negative  as 
well  as  positive.  We  were  able  in  this  way  to  obtain  the  fundamental 
equations  of  Luttinger  and  Ward  in  a  new  way,  in  particular  LW  '  s 

expression  for  the  grand  partition  function  is  obtained  in  a  concise  way. 

* 

The  analogy  between  time  and  temperature  Green's  functions  while  close 
is  not  perfect.  For  example,  in  the  energy -momentum  representation 
(henceforth:  Ep-rep)  the  analogy  breaks  down  inasmuch  as  the  retarded 
and  advanced  temperature  Green's  functions  in  this  representation  are 
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identical,  whereas  this  is  not  true  for  the  time  Green's  function.  At 
this  point  it  is  shown  (Chapter  IV)  that  there  is  a  natural  way  to  split 
the  temperature  Green' s  function  into  two  parts  which  are  separately 
related  to  the  Ep-reps  of  the  advanced  and  retarded  time  Green's 
f unc tions  . 
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CHAPTER  II:  MATHEMATICAL  PRELIMINARIES 

(a)  Introduction. 

In  this  chapter  we  present  a  brief  outline  of  that  aspect  of 
the  formal  theory  of  Green's  functions  which  pertains  to  the 
expositions  and  methods  of  field  theory.  In  section  (b)  the  formalism 
of  the  functional  diffentiation  is  outlined;  we  will  use  this  formalism 
to  define  the  Green's  functions  appropriate  to  our  proolem. 

Section  (c)  contains  a  brief  resume  of  Schwinger's  variational 
principle,  which  is  used  to  derive  an  explicit  expression  for  the 
Green's  function.  In  section  (d)  the  salient  results  are  summarized  in 
a  form  that  will  be  convenient  in  the  succeeding  chapters. 

It  might  be  appropriate  at  this  point  to  make  some  descriptive 
remarks  about  the  basic  mathematical  problem  with  which  we  are 
concerned  in  this  chapter. 

In  the  ordinary  context  of  Green's  functions  one  encounters  linear 
differential  operators  and  sources;  the  latter  are  completely  specified 
at  the  outset.  In  the  theory  of  interacting  particles,  however,  one  is 
faced  with  integro -differential  equations  (without  external  sources) 
which  can  be  regarded  as  non-linear,  homogeneous  equations  for  the 
Heisenberg  operators.  Naturally,  there  are  serious  difficulties  in 
extending  the  Green's  function  approach  to  this  pro  lem.  The  tecnnique 
which  will  be  elaborated  in  this  chapter  (Kato  et  al.  ,  19^0)  is  to 
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generalize  the  definition  of  the  Green's  function  for  linear  operaiors 

to  this  non  -linear  case.  In  essence,  the  me.  hod  is  to  introduce  an 

auxiliary  source  as  matter  of  convenience  -  at  the  end  of  the 

A 

calculation  a  limit  is  taken  in  which  this  source  is  set  equal  to  zero. 
The  generalized  Green's  functions  thus  obtained  provide  a  convenient 
vehicle  for  expressing  various  quantities  of  physical  interest. 


(b)  Functional  Differentiation 

The  formalism  of  functional  differentiation,  which  has  been 
discussed  by  many  authors  (e.g.  Schweber,  19  1  pp.  193;  Katoet.  al.  , 
19  0;  Akhiezer  and  Berestezki,  1953  pp.  10),  provides  an  elegam 
formal  sm  for  the  expression  of  formulae  associated  with  tne  Green  s 
function  technique.  This  formalism  provides  a  transparent  "natural" 
general  zation  of  the  Green's  function  technique  to  the  methods  of  field 

theory . 


respect 


The  functional  derivative  o  he 
co  (p(X)  -  denoted  by 


expression 


7  [  (j)(x) j  wi,-h 

-  is  defined  by 


{ s?[0<*)]  ,  .  =  iim 

J  d<pd)  r  s 

where  j-  (x)  is  reas onably  arbitrary  (e.g.  a  "s moot. n"  and 
integrable  function  of  x) ,  and  we  can  assume  the  domain  of  integra  ion 
spans  'he  region  of  definition  of  [_$(*)]  *  i0H0Vv  A 

pertinent  formulae  follow  immediately  from  tne  definition  1). 
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&(p(x)  r  .  / 

— =  £(*“*) 
(5  </>(*) 


S  J  QU)  4(x)dx 

8  cp(x) 


=  QM 


(2a) 


(2b) 


The  Green's  function  associated  with  the  equation 


L  (f(x)  =  J  (x)  )  (3a) 

plays  the  role  of  the  inverse  of  L  in  the  sense 


In  equations  (3)  it  has  been  assumed  (and  we  shall  continue  to  do  so 
throughout) that  L  is  non -singular  (i.e.  there  are  no  non -trivial 
solutions  to  the  equation  L,  (j)  —  0  consistent  with  the  boundary 
conditions).  Homogeneous  boundary  conditions  at  infinity  are  assumed 
throughout  and  the  integration  in  (3b)  is  to  be  carried  out  over  the  entire 
3 -dimensional  space.  Equation  (3b)  serves  as  the  de:inition  oi  the 
Green's  function  associated  with  equation  (  ^a)  which  satisfies  me  above- 
mentioned  boundary  conditions. 


One  can  readily  convince  oneself  that  for  a  linear  operator  in  (  _a) 
the  solution  (3b)  is  such  that  the  integral  in  (3b),  depends  linearly  on 
the  source  function  J(x)  -  or  what  amounts  to  the  same  thing,  G  x,  x)  is 
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is  independent  of  J(x)  for  linear  operators.  This  observation  along 
with  formula  (2b)  leads  to  an  alternative  definition  for  the  Green's 
function  for  a  linear,  non-singular  operator: 


£(*>■*)  - 


Sltf) 


(4) 


In  order  to  generalize  this  formula  to  the  case  of  non-linear  operators 
we  make  use  of  the  condition  that  the  operators  are  non -sing ular ;  i.e. 
that  (j l  0  as  J  0  .It  can  be  shown  that  in  the  non-linear 

case  the  solution  (p  can  be  expanded  in  a  "Taylor -like "  series 
about  T  =  0  ,  as  follows  (Kato  et.  al.  ,  I960): 


<t>w  = 


r 


j 


h 

l  Ji*) 


jU)d* 

J  =° 


r 


+ 


i 


iM 


j 


J{x)  W 


where  we  have  used  (j)(X )  —  0  .  That  this,  in  fact,  is  a  "natural" 

7=o 

generalization  of  the  familiar  Taylor  series  can  be  seen,  for  example, 

by  considering  the  first  term  in  the  above  series  in  the  case  where 

usual  partial  differentiation  is  permissible.  Thus,  if  we  regard 

0(pc)  ,  which  is  actually  to  be  regarded  as  a  functional  of  J,  as  a 

known  function  of  J,  the  expression  (5)  is  just  the  expansion  of 

0  (  J  )  about  the  "point"  J  =  0,  and  the  first  term  corresponds  to 

ld(j>/\-)  ’  J  (x)  as  one  would  expect  for  a  Taylor  series.  The 

V  t  d  J /j_0 

expansion  (5)  is  useful  if  J(x)  is  small  enough  to  insure  convergence.  In 
the  case  where  the  source  J(x)  is  infinitearrally  small  it  is  expected  that 
only  the  retention  of  the  first  term  in  (5)  is  necessary  (Gottfried  and 
Pieman,  I960).  In  this  case  we  recover,  for  the  non-singular  and 


13 


non-linear  operator,  equation  (4)  in  the  following  form: 


=  iM 

'  (6) 

J  =0 

We  remark  in  summary  that,  for  non-linear  operators  (e.g.  for 
all  non-trivial  interaction  operators  in  field  theory)  which  are  also 
non-singular,  equation  (6)  provides  a  relatively  simple  formula  for 
the  Green's  function  for  weak  (or  at  least  for  infinitesimal)  sources. 
One  should  note  at  this  point  that  in  the  physical  problem  to  be 
considered  here  the  "source"  should  be  regarded  as  an  external 
potential . 

The  transcription  of  the  Green' s  function  formalism  to  field 
theory  is  most  conveniently  done  when  the  latter  is  considered  in  the 
so-called  Heisenberg  picture,  wherein  the  dynamics  of  the  physical 
quantities  are  contained  in  the  operators  themselves.  The  physical 
reasons  for  this  are  discussed  in  Schweber  (1961,  pp.  651);  for  our 
purposes  we  merely  remark  that  in  field  theory  the  analogue  to  the 
condition  (j>  0  as  J  0  ,  as  will  be  shown  below  (pp.lf,lG),  is 

satisfied  in  the  Heisenberg  picture  in  a  direct  way.  Now  in  the 
He  1  senberg  picture  the  equation  of  motion  of  the  absorption  operator 
lj/ (x)  ,  i  .e.  of  the  operator  that  removes  a  particle  at  the  point 

x( x  =  r  ,  t  ),  is  given  by 


1 
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f 

(cBt  +  \^+(x){J[xix.)^'(x)d)c  l^ix) 


(7) 


wherein  U(x  x)  is  the  potential  of  the  particle  pair  interaction  and 
we  have  set  f)  -  £  m  =  1  .  This  equation  is  a  homogenous  (non¬ 

linear  )  equation  for  \jy  (x)  and  should  be  compared  with  (3a)  with 
J(x)  =0.*  It  is  convenient  and  customary  (Kato  et.  al.  ,  I960)  to 
replace  our  interest  in  ^(x),  in  equation  (7),  by  a  consideration  of 
•  i-e-  t^ie  exPectati°n  values  of  "^(x). 


We  remark  here  that  this  replacement  has  "physical  sense" 
inasmuch  as  only  the  expectation  values  of  an  operator  are  meaningful 
physically.  However,  in  the  presence  of  a  source  we  are  still  faced 
with  the  difficulty  of  defining  / j  uniquely,  since  the  introduction 
of  a  source  into  the  Lagrangian  (Schwinger,  19  51)  has  the  result  that 
the  state  vector  of  the  system  (the  eigenvector  of  the  Hamiltonian)  is 
no  longer  stationary  (we  are  interested  in  a  nonstationary  source -an 
external  perturbing  potential).  This  difficulty  is  overcome  by 
considering  the  source  to  be  time  independent  (for  convenience  we 
take  it  to  be  zero)  at  the  beginning  and  at  the  end  of  the  process  under 
scrutiny.  With  these  conditions  on  the  source  in  mind,  the 
expectation  value  of  1 j/  can  be  defined  in  a  more  general  way  than 
indicated  above. 


*  The  bar  on  l^(x)  is  used  as  a  reminder  that  we  are  dealing 
with  Heisenberg  field  operators. 


15. 


In  fact,  we  will  consider  the  following  as  the  fundamental  definition 
for  the  expectation  value  for  any  Heisenberg  operator  Q(x): 


<Q>  =  Tr.{pQ}/Tr.{p]  , 

where 

p  *  \™f>  £*P[-P  £m  +  °<  Nj  <Cmi  I  , 

m 

m-  being  an  initial,  and  a  final  state  vector  of  the  system  in  the 

Heisenberg  picture.  A  detailed  discussion  of  the  meaning  of  the 
summation  is  given  in  Appendix  A.  In  the  above  definition 
(equations  (8)  )  we  have  specialized  to  a  system  in  a  grand  canonical 
ensemble  which  is  the  only  one  considered  in  this  dissertation.  In 
equations  (8),  E  ^  is  the  energy  eigenvalue  of  the  m-th  state  with  N 
particles  present,  (3  '  =  kT,  k  being  Boltzmann 1  s  constant,  and 

o<  is  related  to  the  chemical  potential  through  =  °/4  .  Note 
that  if  the  spectrum  of  initial  and  final  states  is  identical,  then  jD 
is  the  usual  trace  of  the  density  operator  for  a  grand  canonical  ensemble 
(see,  e.g.,  ter-Haar,  19 5 pp.  149). 


(8a) 


(8b) 


We  now  focus  our  attention  on  (spin  indices  are  suppressed 


throughout) 


b  < 


J=0 


<5 


(9> 


' 
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This  equation  defines  the  Green's  function  of  interest  to  us  in  this 


0 


as  J 


investigation.  We  note  immediately  that 


i.e  the  non-singular  condition  is  satisfied  because  of  the 
properties  of  the  absorption  operator  \j/  (x)  and  the  orthogonality  of 
the  Heisenberg  state  vectors.  It  is  of  interest  to  mention  that  the 
discussion  presented  here  (Kato  et.  al.  ,  I960)  is  also  of  current 
theoretical  interest  in  problems  involving  the  so-called  "response" 
of  a  system  to  infinitesimal  disturbances,  e.g.  in  the  study  of  the 
stability  of  physical  systems  (Gottfried  and  Pieman,  I960). 

In  order  to  get  a  manageable  expression  for  the  left-hand  side 
of  equation  (9),  appeal  must  be  made  to  the  physical  laws  governing 
the  variation  of  Heisenberg  operators.  A  convenient  variational 
principle  for  this  purpose  is  Schwinger’s  dynamical  principle  which 


is  outlined  in  the  next  section. 
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(c)  Schwinger's  Dynamical  Principle. 

The  dynamical  principle  of  Schwinger  (19  51,  1953)  is 
summarized  by  the  equation 


r 


0< 


=  iU\b  XcTx 
X  J 


(10) 


where  is  the  Lagrangian  density  of  the  problem,  j  oC^>  is  the 

final  and  J  the  initial  state  vector  of  the  system.  In  words, 

equation  (10)  states  that  a  small  variation  of  the  external  sources 
(i.e.  of  appropriate  terms  in  X  )  results  in  a  variation  of  the 
transition  probability  amplitude  <^o( 


The  prescription  (10)  was  shown  by  Schwinger  (1951,  1953) 
to  be  the  generic  formula  for  all  the  fundamental  equations  and 
relations  in  quantum  mechanics.  In  particular,  a  simple  demonstration 
of  the  equivalence  of  this  prescription  to  the  Feynman  formulation  of 
quantum -mechanics  (Feynman,  1948)  is  outlined  in  the  review  paper 
of  Kato  et.  al.  (I960). 


Note  that  the  transition  amplitudes  <^c<  j  in  equation  (10) 
refer  to  pure  states  and  are  distinct  from  the  canonical  averages  of 
Heisenberg  operators  which  occurred  in  equation  (8).  Here  |  0( 
and  |  are  eigenvectors  of  the  Hamiltonian.  In  the  present 

investigation  we  are  concerned  with  a  system  which  is  characterized 
by  a  given  temperature  and  whose  particle  number  can  vary.  Thus 


we  cannot  character ;  ze  the  system  by  an  eigenvector  of  the  Hamiltonian 
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either  before  the  introduction  of  the  source  or  after  the  source  is 
switched  off.  In  fact  this  was  the  reason  that  we  used  statistical 
expectation  values  as  defined  by  equations  (8). 

Naturally , a  'source  term'  in  an  equation  for  ^  must  be 
considered  as  deriveable  from  a  change  in  the  Lagrangian  density 
of  the  Fermi  field.  Thus,  in  addition  to  the  Lagrangian  appropriate 
to  the  actual  system,  we  include  a  term  of  the  form 

where  J(x)  is  an  arbitrary  external  perturbation  which  anticommutes 
with  j^(x).  With  this  change  in  the  Lagrangian  density  we  get  from 
equation  (10),  the  following  equation  (the  proof  and  references  are 
given  in  Appendix  A): 


G(x>x) 


where  the  angular  brackets  are  to  be  interpreted  by  equations  (8),  and 
where  T  is  Wick's  chronological  operator  (Wick,  19  50): 


for  t>t‘ 

—  —  W  (x  J 

for  t  <  t  , 

In  the  next  section  the  results  to  this  point  are  summarized  and 
an  expression  for  the  Gibbs  potential  in  terms  of  the  Green  s  function 


i  s  written  down. 
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(d)  The  Equation  of  Motion  for  the  Green's  Function. 

We  start  with  a  brief  summary  of  the  results  obtained  so 
far.  The  usual  Green's  function  for  linear,  non-singular  differential 
operators  which  act  on  a  scalar  function  (j)  (x)  (i.e.  (p  (x) 

is  a  one  component  function  which  is  scalar  with  respect  to  rotations 
in  ordinary  3-dimensional  space)  in  the  presence  of  a  source  J(x), 
was  given  (Kato  et.  al  ,  I960)  in  terms  of  functional  differentiation: 


S  (p(*> 
h  Jib 


( Linear  case) .  (4) 


The  case  of  non-linear,  non-singular  operators  operating  on  a  scalar 
function  with  infinitesimal  scalar  sources  was  shown  to  lead  to  an 
analogous  equation: 


G(X>X)  - 


S  0 (■*) 

$  J  I*) 

J  =0 

(Non-linear  case).  (6) 


In  the  derivation  of  formula  (6)  it  was  assumed  that  the  solution  } 

<P(x)  )  depends  linearly  on  the  source  (we  have  simply  neglected 
the  higher  order -terms  --  which  is  presumably  justifiable  for 
infinitesimal  sources).  At  this  point,  rather  than  considering  the 
Heisenberg  operator  (x),  we  consider  the  equation  for  , 

where  the  expectation  value  was  defined  to  be  the  grand  canonical 


average  (equation  8).  (It  should  be  mentioned  here  that  a  formulation 


20 


based  on  petit  ensemble  average  has  been  published  'Brout  and 
Englert,  19f0)  --  however,  'he  petit  ensemble  approach  to  the  problem 
is  much  more  complicated  than  that  of  the  grand  ensemble.  The 
Green's  function  was  defined  to  be 


-  (j  (■*>*)  = 


^  <  b-  (A!)> 

5j<*>  l 


(9) 


which  is  a  "natural"  generalization  of  equation  (fc)  ,  Schwinger  s 
dynamical  principle  was  then  utilized  to  obtain  an  explicit  formula 
for  Gtx,  x)  in  terms  of  the  Heisenberg  absorption  and  creation 
opera  ors  and  1/  (X) 


G(x>*)  =  i  y  •  an 

To  obtain  the  equation  of  motion  sat  sfied  by  the  Green  s 
function  we  now  proceed  as  follows  Consider  the  so-called 
Schwinger  identity  (Klein  and  Prange,  1958): 

dt  (X)j  =  I  jdt  :f  lx)f  +  k  (X-X).  (12) 

(This  equation  is  proven  in  a  simple  way  in  Appendix  A}  Equations 
(12)  and  (11)  coupled  with  the  equations  of  motion  of  the  Heisenberg 
operator , 

=\yto,H]_  ,  (U) 
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lead  to  the  equation  of  motion  for  the  Green's  function  (H  is  given 
by  equation  21) : 


+  Vc‘)  fr(x.x) 


(14) 


The  expression  in  angular  brackets  can  be  simplified  in  the  following 
manner.  We  define  a  dimensionless  parameter  A  (e.g.  see  Klein 
and  Prange,  1958)  which  characterizes  the  interaction,  (J  (  Xix) 

A  U[X>X )  ■  Furthermore  we  define  the  "irreducible  self 

energy  operator"  (Klein  and  Prange,  1958;  Klein,  I960;  Klein,  1961) 
by  the  equation 


r 


r 


cCx  M (*>*)  &(*,*)  - 1 A 


A  ^ (h y (*) j/ tr)jf  f(x)\  d  •  ( 1 5) 


Since  our  discussion  will  eventually  lead  to  a  perturbation 
expansion,  it  is  useful  to  define  the  "non-interacting"  Green's 
function,  (X,X)  .  This  function  is  defined  by  the  equation: 

(i  3t  +  Vr  )  (jC>U>x )  =-&  {*-*)  »  (16) 

i.e.,  the  equation  obtained  by  setting  (J  —  0  in  (14).  We  now 
insert  (15)  into  (14)  and  rewrite  the  resultant  equation  in  the  form 
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The  definition  (16)  can  now  be  used  to  convert  (17)  into  an  integral 
equation.  The  result  is  (Klein  and  Prange,  1958) 


G(x,i) 


(18a) 


This  equation  can  be  written  in  the  more  suggestive  form 


G-  -  GM  +  G  i0)M  G 


(18b) 


(18c) 


The  last  expression,  (18c),  is  seen  to  be  a  perturbation  expansion; 
in  fact,  as  we  shall  see  in  Chapter  III,  it  bears  a  close  similarity  to 
equation  (38)  in  the  Luttinger  and  Ward  paper. 

We  pause  here  for  a  general  qualitative  remark  on  the  importance 
of  the  assumption  that  our  problem  can  be  solved  by  perturbation 
procedures  (i.e.,  by  terminating  the  expansion  (18c)).  The  physical 
significance  of  this  assumption  can  be  seen  in  that  it  rules  automatically 
against  the  possibility  of  a  bound  fermion  pair  of  the  Cooper  type  (Cooper, 
1956).  This  is  because  the  energy  gap  term  (  in  Cooper's  notation) 
can  be  obtained  from  the  methods  of  field  theory  only  by  summing  an 
infinite  number  of  terms  in  an  equation  like  (18c).  Thus  the  summation 
indicated  by(l8c)can  be  written  formally  as: 

/ 

G  =  — - -  , 

[GkT  -  M 


wherein,  crudely  speaking. 


m 


Cooper  s  theor 


Y  and  [Q 


M  plays  the  role  of 

! 

is  the  operator  inverse  to  the 


non-interacting  Green’s  function.  At  the  same  time  it  is  well-known 
that  the  accepted  theory  for  superconductivity,  the  so-called  BCS 
theory  (Bardeen  et.  al  ,  1957)  rests  on  the  existence  of  "Cooper 
states  ",  Thus,  in  this  case,  the  usual  perturbation  procedure  is 
unable  to  unravel  the  physical  situation.  In  general  it  is  to  be 
expected  that  the  energy  disconrinuities  and  other  processes  involving 
phase  changes  (such  as  was  suggested  by  Tramor  and  Dixon  for 
nuclear  matter,  (1955))  cannot  be  explained  by  a  perturbation 
expansion.  The  fact  that  perturbation  techniques  are,  nonetheless, 
widely  used  in  many  body  theory  expresses  in  part  the  difficulty  of 
those  problems  from  both  mathematical  and  conceptual  points  of  view; 
on  the  other  hand,  it  is  relatively  easy  to  assign  a  pictorial  meaning 
to  the  perturbation  terms  involved.  Apart  from  this,  perturbation 
expansions  are  useful  because  the  very  nature  of  their  failure  in  a 
given  problem  helps  to  delineate  those  features  which  are 
characteristically  important  for  that  problem.  As  an  example 
consider  once  again  the  theory  of  superconductivity:  if  one  attempts 

a  perturbation  expansion  for  the  electron  interactions  in  metals,  the 
second  order  terms  have  a  singularity  at  the  fermi  surface.  The 
failure  of  perturbation  techniques  led  eventually  to  Cooper's  theory 
which  provided  a  sound  basis  for  the  description  of  the  phenomenon  of 
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The  expression  for  the  energy  of  a  system  in  a  grand  canonical 
ensemble  is  given  by  (cf.  ,  e.g.,  ter  Haar,  1955  pp.  149  ) 


c  _  Jh 

where  the  last  identity  expresses  our  general  definition  for 
expectation  values  (equation  8).  (Note  that  this  quantity  is  evaluated 
with  an  infinitesimal  source  J  and  then  the  source  is  set  equal  to  zero.) 
The  chemical  potential,  M  ;  is  defined  as  usual  by  the  requirement 
that 

j  /  he  average  number  of  particles;  -  <Vv  J>  •  <2°) 


The  Hamiltonian,  H,  in  the  Heisenberg  picture  can  be  written  as 

(J  *  0) 


wherein  terms  enclosed  by  colons  are  to  be  taken  in  the  "normal"  order 
(Wick,  1950)  ;  i.  e.  ,  in  any  sequence  of  operators,  all  referring  to  the 
same  time,  absorption  operators  are  to  be  brought  to  the  right  of 
creation  operators.  (We  assume  throughout  that  the  interaction  is 
instantaneous  :  i .  e  .  ,  two  particles  interact  only  at  the  same  time  .  ) 
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We  shall  not  discuss  the  reason  for  having  the  operators  in  their 
normal  order  (Wick,  1950),  but  merely  remark  that  in  this  way 
some  unobservable  divergences  are  automatically  removed  from  the 
theory  at  the  outset  (Dyson,  1952,  pp.  156). 

In  closing  this  review  chapter,  we  would  like  to  outline  the 
derivation  of  an  alternative  to  the  usual  expression  for  the  Gibbs 
potential.  The  Gibbs  potential,  A,  is  usually  defined  by 


Tr.  {-&rp[-/3 (H-/< fi/jjj  =  -£*7>[/Mj  ,  (22a) 


or 


A  =  L 


(22b) 


The  alternative  expression  for  A  which  is  of  interest  to  us  is 
obtained  as  follows:  we  characterize  the  interaction  by 
the  replacement  .  U(x,x)  A  U(x,  x)  and  thus  write 


H  ~  Ho  +  'X  Hi  ) 


« - 1 


( 


)\  ,  we  obtain  (cf. 

9/1  ' 


and  upon  differentiation  with  respect  to 
Klein,  1961a) 


dP 


H, 


p>l/ 


(23) 
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Converting  (23)  into  an  integral  equation,  it  reads 


(24) 


This  equation  leads,  via  the  definition  equation  (15),  to  our 
alternative  expression  for  the  Gibbs  potential: 


This  equation  was  obtained  by  Klein  (1961b). 

It  is  important  to  note  that  A  /  involves  the  chemical 

1  I  >  =  0 

potential  as  determined  by  equation  (20).  That  yK  ,  the  chemical 
potential,  might  be  different  from  the  chemical  potential  of  free 
fermions  (labelled  by  M0  )  forms  a  cornerstone  in  the  LW  scheme 
for  elimination  of  the  so-called  anomalous  diagrams. 


In  the  next  chapter  it  w:ll  be  shown  that  equation  (25)  leads  to 
the  basic  equation  (45)  in  LW  '  s  paper  for  the  Gibbs  potential, 
although  (25)  was  derived  by  considering  temporal  developments  while 
LW  used  the  BDD  approach  (cf.  p.  3). 
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CHAPTER  III:  TIME  AND  TEMPERATURE 

GREEN'S  FUNCTIONS 

Intr  oduction . 

In  the  previous  chapter  we  outlined  the  formal  derivation 
of  the  expression  for  the  single  particle  Green’ s  function.  We 
shall  refer  to  this  expression  (equation  (9))  as  the  "time" Green '  s 
function.  It  is  well  known  that  (c  f .  ,  e.g.(  Thouless,  1961)  the 
"time" Green '  s  function  lends  itself  to  a  relatively  simple 
interpretation  if  free  particles  are  considered  and  complications 
due  to  statistical  averaging  are  deleted.  In  this  case  one  can 
consider  the  non-interacting  "time,,Green '  s  function  as  the 
scalar  product  of  the  ground  state  with  the  state  which  is  formed 
from  the  ground  state  by  the  addition  of  a  particle  a*  the  time  t  , 
allowing  the  particle  to  propagate  freely  until  t'  ,  and  then 
removing  a  particle  from  the  system.  (In  the  above  interpretation 
we  have  assumed  that  the  creation  operator  was  to  the  right  of  the 
absorption  operator;  it  is  clear8  however,  that  a  similar  inter¬ 
pretation  can  be  presented  for  the  reverse  situation.)  It  is 
pos  sible  to  show  that  ( c  f .  s  e  .  g .  ,  Thoule s s  loc  .  cit . )  the 
expectation  value  of  any  single  particle  operator  can  be  obtained 
directly  from  the  "time^Green '  s  function.  We  shall  not  attempt 
the  interpretation  for  the  interacting  Green  s  function  in  a 
statistical  ensemble;  however,  this  can  be  done  by  "stretching 
somewhat  the  above  discussion. 
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We  are  unable  to  give  an  interpretation  for  the  "temperature" 
Green' 8  function  which  is  alluded  to  in  the  title  of  this  chapter. 

Thus  it  must  be  considered  as  a  formal  device  wherein  one 
encounters  what  we  will  refer  to  as  a  "temperature  variable.  " 

This  variable,  which  is  customarily  denoted  by  u,  has  the  dimension 
of  reciprocal  temperature  (i.e.  reciprocal  energy)  and,  in  our 
formalism,  is  bounded  by  -  { 1 ^  U  ^  +  f^/z  where 

ft-'  =  kT. 

The  parallelism  between  the  roles  of  the  temperature  variable 
and  the  time  in  the  Green's  function  formulation  of  the  temperature 
dependent  MBP  was  first  discussed  by  Matzubara  (1955)  who 
considered  systems  with  a  fixed  number  of  particles.  Matzubara 's 
approach  was  extended  by  Fradkin  (1959)  and  AGD  (1959)  to  a 
system  in  a  grand  ensemble. 

On  the  other  hand  it  was  the  per turbation -like  approach  of 
BDD  (1959)  which  was  used  by  LW’  (I960)  for  their  discussion  of  the 
correct  generalization  of  the  Goldstone  (1957)  expansion  for  systems 
at  finite  temperatures. 

Klein  (I960,  1961a)  obtained  the  generalization  of  the  Goldstone 
expansion  for  systems  at  finite  temperatures  by  using  the  "time" 
Green's  function  approach.  Furthermore  Klein  (1961b)  was  able 
to  derive  Landau's  semiphenomological  theory  of  Fermi  liquids 
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(Landau,  19  57)  by  an  extensive  use  of  the  "time"Green 1  s  function. 

In  this  chapter  the  temperature  dependent  MBP  is  reformulated 
in  terms  of  a  new  variable  T  which  will  be  shown  to  play  the  role 
of  either  the  time  (t)  or  the  temperature  variable  (u);  this  new 
variable  replaces  the  time  in  the  theory  outlined  in  the  previous 
Chapter.  In  this  way  the  intimate  interchangeability  between  t  and 
u  is  shown  explicitly.  Further,  our  formulation  permits  the 
assignment  to  u  of  negative  as  well  as  positive  values  and  thus 
enhances  its  similarity  to  the  time.  We  prove  the  validity  of  our 
formulation  by  deriving,  using  the  "temperature11  Green 1  s  function, 
the  fundamental  equations  of  LW  for  the  Gibbs  potential  and  the 
single  particle  propagator. 


(b)  Mathematical  Preliminaries . 

Our  approach  involves  the  following  substitution: 


(1) 


where  t  is  interpreted  as  the  usual  time  variable  while  u  is  to  be 
regarded  as  a  "temperature  variable'1  with  the  dimension  of 
reciprocal  temperature,  Le.  reciprocal  energy.  This  substitution 
has  the  following  consequences:  - 
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(i) 


an 

where  "^/  (r,  t)  is, arbitrary,  reasonably  well  behaved  function, 
for  example,  any  function  which  has  a  Fourier  expansion: 


f(cj) 


1  ? 

a F(r)  e 


Ofia 


I  +  ip* 


r 


(ii)  If  T"  is  subjected  to  a  similarity  transformation  through  a 
2x2  matrix, 
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then  it  follows  immediately  from  formula  (2)  that 


^(r.t)  =  S  S 


(iii) 


(4) 


This  formula  follows  in  an  obvious  way  from  equation  (2).  Its 
consistency  can  be  seen,  for  example,  in 


(iv)  We  now  regard  T  as  a  variable  conjugate  to  the  Hamiltonian; 
hence  for  Heisenberg  operators  such  as  ^  (  f  j  T)  we  have 


> 


(5) 
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where, 

of 


in  addition  to  normal  commutation, 
ijy  and  //  is  implied. 


-  multiplication 


(v)  The  chronological  operator  T  operating  on 
is  now  interpreted  as 


((>) 


wherein  each  matrix  element  the  operation  of  T  is  the  usual  one 
(Wick,  1950);  thus  for  the  field  operators  'U'  and  : 


U,  n  lc,  f') 


1 


and 


ft?,'") 


=  fit  A)  tit 

=  -f(dt)?(r,t)  i<t  , 

]  -  'jf-  tr.*’1*)  u  >  u 

-  -  ^  Vr ,‘U)  “)  u  <  u  • 
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(c)  The  Equation  of  Motion  of  the  Green1  s  Function. 

(i)  The  Schwinger  Identity: 

The  Schwinger  identity  can  be  obtained  in  the  usual  way. 
From  the  definition  of  T,  we  have 


=  H^-T)  jf  (t>T)  O(r-i)  j'  (W)  'jlC'T)  • 


Here 


91  r )  = 


$(t)  o 


0  @(iu) 


whe 


re 


1 

t  >  0 

9  (w  = 

0 

t  <0  > 

i 

H  >  0 

61c  u)  - 

0 

U<  U  . 

Hence  (cf.  Appendix  A) 


dJlf(£T)f(Cj)]=  c>(r-r')S(r-r)  +  /  [drf(r,r)f*(r,T)j  , 

where  we  have  used  the  anti -commutation  relation  for  the  Heisenberg 
operator  \j/  O-ind  j/  > 

(ii)  The  Equation  of  Motion  of  the  Heisenberg  Operator  ^  • 

Using  equation  (5)  we  get 


idry  =  -  py  +  ijd V  U(if-rl):^lrT)^(r  bj/(c,r): 


A 


5  (7) 
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where  H  is  given  by  equation  (21)  of  Chapter  II.  We  note  that  the 
operators  involved  in  (7)  all  refer  to  the  same  T  .  The 
meaning  of  the  colons  is  the  same  as  in  Chapter  II  (equation  21)  ; 
however,  here  all  the  formulae  are  matrix  relationships. 

(iii)  Definition  of  Expectation  Values  in  a  Grand  Canonical  Cnsemble. 

In  the  formalism  used  here  the  Heisenberg  state  vectors 
must  be  considered  as  having  two  components.  Thus,  our  (2x2)  matrix 
operator  has  two  eigenvectors  corresponding  to  each  eigenvector  of 
the  usual  Hamiltonian.  We  shall  use  the  notation 


Here  m  denotes  the  Heisenberg  state  vector:  j-j  j  my  =  e  I  my 
and  /V0  I  m y  -  /  my  where  /-/  and  /\!0 

are  the  usual  Heisenberg  Hamiltonian  and  particle  number  operators. 
Either  one  of  the  above  states  (i.  e.  ,  f  or  |  )  can  be  used  to 

obtain  all  the  properties  of  the  system  as  will  become  clear  below. 
However,  to  retain  the  generality  of  our  approach  we  will  use  the 
notation 


m 


m 


o  \ 


o 


m 
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however  all  results  should  be  interpreted  as  either  f  or  jf 
In  this  way  the  formulae  can  be  developed  in  a  unified  fashion.  W  ith 
this  preamble  we  can  proceed  with  the  pertinent  definitions.  The 
expectation  value  of  an  arbitrary  operator  Q  is  defined  by  its 
grand  canonical  average:  - 


(8a) 


with 


(8b) 


and 


CO 


m 


-  ji  E  +  o(  N 


J 


(8c) 


where  the  terms  are  defined  in  the  same  manner  as  in  conjunction 
with  equation  (8)  of  Chapter  II. 


(iv)  Equation  of  Motion  of  the  Single  Particle  Green's  Function: 

The  single  particle  Green's  function  can  now  be  defined 
in  the  usual  way:  - 


t 


&(x>x) 


L 


(9) 
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where  X  —  ^  C  >  L  )  * 

Using  equations  (6)  (7)  (8)  (9)  we  get(cf.  the  derivation 

of  equation  (14)  in  Chapter  II)  the  equation  of  motion  of  the  Green's 
function: 

(10) 

(11) 


=  b(x-x)  -  ‘  K  (  x,x) 

where 


which  is,  of  course,  formally  equivalent  to  equation  (14)  of  Klein 

and  Prange  and  equation  (14)  of  Chapter  II.  The  equation  of  motion 

•  /  /  \ 

of  the  non-interacting  single  particle  Green's  function,  (j-  (  *>*)  , 

is  given  by 

-<  3r  -  =  &(*-*)=  S(t-t) £(r-r)  ,  az) 


i.e.  by  setting  X  —  0  eTuatl°n  (I1  )• 

(v)  The  Iterative  Equation  for  the  Green's  Function. 

As  in  Chapter  II  (section  (d)  we  define  the  "irreducible 
self  energy  operator"  by  the  equation: 
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... 


i'x  ^T|d3r  : 


(13) 


Substituting  (13)  into  (10),  the  latter  can  be  converted  into  an 
integral  equation  upon  utilization  of  (12).  The  result  is 


r 


G  ( *.i)  -  Gb)  (*.i)  +  I ,  04.) 


in  the  more  suggestive  form 


As  before  we  can  write  equation  (1-ia.j 


G 


=  G 


{0) 


(14b) 


Equation  (14b)  is  the  iterative  equation  for  the  single  particle  Green’ s 
function. 

In  the  next  section  we  will  prove  the  validity  of  the  equations 
given  in  this  section. 

(d)  Proof  of  the  Validity  of  the  Present  Formulation. 

If  one  regards  the  expectation  values  that  were  considered  in 
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the  previous  section  as  taken  with  respect  to  the  states  designated  by 
j  (i.e.  ;  )  then,  clearly,  all  the  equations  considered  in 

that  section  become  identical  with  those  considered  in  the  previous 
chapter.  Hence,  this  case  requires  no  further  discussion.  The 
validity  of  our  equations  does  require  special  attention  when  we  consider 
the  expectation  values  as  taken  with  respect  to  the  states  designated 
by  |.  (i.e.  )  of  the  previous  section.  The  verification  of 

the  expressions  for  this  case  will  be  provided  below  by  utilizing  our 
equations  to  derive  the  fundamental  equations  of  LW  for  the  Gibbs 
potential.  Having  established  this  equivalence,  one  can  then  turn  to 
the  analysis  of  Luttinger  and  Ward  for  a  physical  interpretation  of 
the  mass  operator. 


(i)  The  Temperature  Green's  Function  for  a  Free  Particle  in  the 

Energy -Momentum  Representation. 

The  free  particle  Green's  function  is  obtained  by  setting 
G  -  0  in  the  expression  for  the  interacting  Green's  function,  i.e. 
equation  (12)  Chapter  II,  The  result  is 


It  might  be  emphasized  that  in  the  present  context  :x  (rt  (  u') 

We  now  expand  the  operators  f  (x)  and  ^  <*)  by  Fourier 

series  -  it  is  implicitly  assumed  throughout  that  this  i3  possi  To 
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(cf.  equation  (2),  p.  3  0  ).  We  limit  ourselves  to  consideration  of 

Green's  functions  which  depend  upon  the  spatial  differences 


f  -  G  )  only;  this  is  equivalent  to  the  assumption  of  a 


uniform  medium  (the  significance  of  this  assumption  is  discussed 
below).  We  then  obtain 


r(o) ,  / 

( 7  (*•*) 


where  V  is  the  volume  of  the  system  and  periodic  boundary  conditions 
have  been  assumed.  The  details  of  the  calculation  leading  to  equations 
(16)  will  be  given  in  Chapter  IV. 

We  would  like  to  pause  at  this  juncture  to  discuss  the  assumption 
of  uniform  density  that  was  used  in  obtaining  equations  (16).  The 
possibility  that  self -interacting  systems  of  identical  particles  may 
possess  a  non-uniform  density  has  been  considered  in  the  literature. 
For  example,,  Trainorf  Bhatia  and  Sood  (1958)  considered  a  periodic 
model  for  nuclear  matter  and  were  able  to  calculate  values  of  the 
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nuclear  density  and  the  binding  energy  per  nucleon  in  good  agreement 
with  experimental  values.  The  possibility  of  obtaining  density 
fluctuations  in  many  body  fermion  systems  has  been  recently 
investigated  by  several  authors  (Overhauser  (I960),  Kohn  and  Nettel 
(I960),  Gross  (1960),  Brout  (I960),  Henley  and  Ruijrok  (I960),  Werner 
( 19 61) ,  Thouless  (1961)  who  were  confined, because  of  computational 
difficultieSjto  the  limit  of  weak  attractive  interactions.  However, 
attention  has  been  restricted,  for  the  most  part,  to  uniform  media; 
this  is  essentially  because  of  the  appreciable  difficulties  involved  in 
describing  the  development  of  density  fluctuations  from  uniformity  as 
the  interaction  between  particles  is  turned  on. 

(ii)  Derivation  of  the  LW  Representation  for  G  [v)  (x,x). 

We  are  now  in  a  position  to  obtain  the  propagator  of  LW  from 
our  equations  (16)  using,  in  fact,  a  technique  devised  by  LW  themselves. 
To  accomplish  this,  we  need  to  restrict  the  variable  u  ;  we  impose 
on  u  the  condition  jLlJ  ^  /3/?  where  6  is  the  usual 

inverse  temperature  factor.  Note  that  u  can  take  on  negative  as  well 
as  positive  values.  We  note  that  in  equations  (16),  the  factors  in  the 
terms  which  contain  explicitly  are  analytic  functions  of 

.  it  is  expedient  now  to  analytically  continue  these  factors 
into  the  complex  S?  plane,  which  we  label  by  J  .We  can 
now  represent  the  terms  which  contain  So  explicitly  in  equations 
(16)  as  follows: 
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4Xp  [+  £p(U-U)~j 

/  +  Mf>[p  (tr'S*)] 


27Ti 


&/>[-!  (U-uj] 


5-£r 


4XP 

l+4Xp[-/i  (tfy*)] 


l 

JTi 


J 


expj-j(u-uij 

Hxyp[-/s  O/'jj 


for  U<(I  ,  (|7b) 


where  fl  is  defined  in  Fig.  1: 


plane 

A 

L - - J 

— t - 

- - 

V  -r 

1  d 

Fig.  1:  The  Contour  Defining  f\ 

In  Fig.  1  the  path  is  chosen  to  be  below  (above)  the  poles  of  the 
integrands  in  the  upper  (lower)  half  of  the  complex  ^  plane.  With 
this  choice  of  H0  the  proof  of  these  integral  representations 


follows  trivially  from  Cauchy's  theorem,  there  being  no  poles  on  the 
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real  axis  except  at  .  Now  the  integrands  in  equations 

(17)  are  zero  at  A  and  B  in  Fig.  1  (i.e  at  t.  °°  );  further,  if  we 

analytically  continue  the  integrands,  as  functions  of  f  into 

the  whole  complex  plane  we  note  that  they  are  zero  on  the  large 
semicircles  above  and  below  the  real  axis.  Hence  we  can  replace  the 

V 

contour  P0  by  the  contours  I  which  are  shown  in  Fig.  2 
(This  technique  is  due  to  LW) . 

-f  pi  An  e 

/ 

/ 

'  r 


Fig.  2:  The  Contours  Defining 

In  other  words  we  have  replaced  the  residues  of  the  integrands  in 
equation  (27)  on  the  real  axis  by  their  residues  in  the  whole  complex 
plane  -  except  for  the  real  axis.  Here  the  advantage  of  this 
procedure  becomes  apparent  -  the  equations  (16)  take  on  an  identical 


form: 
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JMp[i?-( C'f)  - 


f, 


(18) 


zAll  jri  +/•  \  /=  o>ti±z.-“ 
ft  x 


With  this  expression  we  can  label  the  single  particle  non-interacting 
Green' s  function  in  the  energy -momentum  representation  by 


(19) 


i.e.  we  consider  the  set 


,  whe  r  e 


0,  (r.’d  s  ^pL;P-r-s,v_ 


(20) 


to  be  a  complete  set  in  the  sense  that  G(x,x),  G  (x,x)  and  M(xsx) 
can  be  expanded  in  terms  of  it  in  the  4  dimensional  range:  Z (S  V 
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This  point  is  discussed  further  in  Appendix  B. 

(iii)  Derivation  of  LW '  s  Fundamental  Equations. 

The  formula  connecting  the  non-interacting  Green's 
function  with  the  interacting  Green' s  function  and  the  "irreducible 
self  energy  operator"  was  obtained  in  section  (d)  (v)  of  this  Chapter: 


&  (*>*) 


J 


(14a) 


Expanding  the  operators  in  this  equation  in  the  set  introduced  above 
(equation  20)  leads  to  (after  the  space  integrations,  which  are  done 
trivially) : 
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Now  the  integrations  of  the  variables  Lit  and  U2 

over  the  range  -  /3/2  to  +  ft/z  can  be  readily  carried  out 

with  the  result: 


With  this,  equation  (21)  takes  the  following  form: 


/ 


and,  because  of  the  linear  independence  of  the 
(deleting  the  indices  £  and  ) 


,  we  get 


£  =  t  > 


(22) 


which  is  identical  in  form  with  equation  (14b),  except  that  here  the 
terms  are  in  the  energy -momentum  representation.  Utilizing 
equation  (19)  we  can  rewrite  equation  (22)  in  form 
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(23) 


Equations  (22)  and  (23)  should  be  compared  with  equations  (38)  and 
(39)  in  the  paper  of  LW. 

We  are  now  in  a  position  to  derive  the  LW  equation  for  the 
Gibbs  potential  (equation  (45)  in  LW's  paper).  Recalling  that  we  had 
(Chapter  II  equation  (25)  )  the  following  expression  for  the  Gibbs 
potential  A, 


(25) 


o 


We  now  expand  once  more  the  functionals  M  and  G  in  terms  of 
equation  (20)  and  get,  after  straight-forward  integration 


(24) 


With  the  aid  of  (19),  (23),  and  (24)  we  are  justified  in 
making  the  following  identifications: 

LW  Notation  Our  Notation 


w 

£/>  (  SjL ) 

The  identification  that  is  sammariued  in  the  table  abov 
can  now  be  used  for  the  physical  interpretation  of  various 
functions  and,  in  particular,  the  so-called  "irreducible  self 
energy  operator."  This  procedure  can  be  done  by  simply 
considering  the  methods  and  definitions  used  by  LW . 

This  completes  our  proof  for  the  validity  oi  our  basic 
equations,  (10)  and  (14),  and  justifies  our  assumptions. 
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CHAPTER  IV:  THE  ENERG Y -MOMENTUM  R EPR ESENTATIONS 

OF  THE  GREEN'S  FUNCTIONS 

Introduction. 

Hitherto  we  observed  the  striking  parallelism  between 
the  time  and  the  temperature  Green's  functions.  This  is  further 
illustrated  below  where  we  show  that  the  definitions  for  the  retarded 
and  advanced  "time"Green' s  functions  find  their  natural 
correspondence  in  the  "temperature11  Green  '  s  functions.  However, 
when  we  consider  the  energy-momentum  representations  (henceforth 
denoted  by  Ep-rep.)  of  these  functions  some  of  the  parallelism  is 
lost. 


The  Ep-reps  for  the  two  cases,  temperature  and  time 
Green's  functions,  are  obtained  in  two  steps.  Firstly,  the  momentum 
representations  are  obtained  by  simply  taking  the  Fourier  transforms 
of  the  3-dimensional  space  dependence  of  the  Heisenberg  operators; 
this  procedure  is  identical  in  the  two  cases.  Secondly,  the  energy 
representations  are  obtained  y  again  by  Fourier -like  transfor mations^ 
with  respect  to  the  time  and  temperature  like  variable  respectively. 
The  resultant  expressions  are  compared  and  discussed.  It  is 
precisely  at  this  point  that  the  parallelism  between  the  time  and 
temperature,  which  was  mentioned  above,  breaks  down.  Now,  it  is 
possible  to  split  the  time  Green' s  in  the  Ep-rep  into  two  parts  (see 
e.g.  Klein  and  Prange,  1958)  in  such  a  way  that  one  part  corresponds 
to  a  retarded,  and  the  other  to  an  advanced,  time  Green's  function. 
The  question  naturally  arises  whether  there  is  any  corresponding 
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division  in  the  temperature  Green's  function  in  the  Ep-rep.  This 
point  is  elaborated  in  the  present  chapter.  In  particular,  the  case 
OO  is  examined  in  detail. 


ft 


(b)  The  Ep-reps  of  the  Temperature  Green's  Function  for 
Finite  Temperatures. 

From  equation  (5)  Chapter  III  we  have 


>  inhere  X  =  (  f  ,  i  li  ]  • 


The  standard  solution  of  this  equation  is 


da) 


Similarly,  one  has 


(lb) 


The  "temperaturellGreen '  s  function  was  given  (equation 


(9),  Chapter  III  )  )  by 
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/T 


Upon  substituting  equations  (1)  into  this  expression,  carrying  out 
a  Fourier  expansion  of  the  spatial  parts,  and,  finally,  demanding 
that  (  ]/r  -f-  Vr'  )  G  =  0  (i.e.  applying  the  condition  for 

uniform  media;  see  comments  of  pp.  3  9,4 0),  we  get 


a  tu)  akiui 

k  k 

-i 


where 

_  a  H  u.  H 

aKu)  =  e  aK  e 


and 


(2) 


(3a) 


(3b) 


and  where  G 


K 


(ueu  )  is  the  Fourier  coefficient  in  the  expansion: 
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G(X’X)  = 


-K 


u^y 


^  ,  /  ck-ic-n  i 

(t„  (u>  U)  6  ^  K 


From  (2)  and  (3)  and  the  cyclic  invariance  of  the  trace  (which  is 
implied  by  the  angular  brackets),  we  have 


=  GK(U~ u)  ~~  6^  ^ 


/ 

tf  '  U'  u,  . 


In  fact,  (2)  can  be  written  more  explicitly: 


y—  -n  fy,a  +  +  v  /  .  ;  -tih  ,  ,  .  v 

G-K(v)  =  -D/fte  \^alaKe  at l"‘<V  - /«-  ^  >  (4 a) 


Aa 


\ 


■  x>Z  e 

fto. 


./S£«>  +c(ll  -  V^A  /  f  lltl 

/M>a\Q^e  CLtlhoSj*  l>o,(4b) 


where 


-/  ^ 
^  =Z  e 

/J<\ 
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and  the  states  j  /v>  Cl  )>  are  defined  by  the  eigenvalue  equation 
H I N \  Ciy  -  j J\}  &y  *  here  N  is  the  number  of  particles 

in  the  state  /J\j)  and  a  stands  for  all  other  quantum  numbers 

that  are  required  to  specify  the  state. 


In  the  previous  chapter  equation  (18)*  and  the  comments  below 
(20)  we  used 


-  oo 


where  we  have  deleted  the  label  £  ^  ^  \LL  I  3  j  LL  \  £  /  7L 

y  -  Lz.A.  .^1 _ jj  Jil  -j-  U  .  From  (5)  we  have  immediately 

X  p  / 


and 


P 

e 


V  >  o  , 


(6  a) 


j 

and 


~  R 

&K  (l] 


-/* 


{or  v<0  . 


(6b) 
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Here  the  superscripts  A,R  signify  the  ’‘advanced*'  (v  >  0)» 

respectively  "retarded"  (v  <  0  ),  temperature  Green's  function. 

In  order  to  utilize  equation$(6)  it  is  expedient  to  rewrite  terms 
occurring  in  equations  (6)  in  a  different  form  (using  a  sum  over 
intermediate  states) 


Thus  equations  (4)  can  be  written  as 


(8a) 


) 


(8b) 


A!na 
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where 


A 


A 

/ 

a i}  a  ft 


5k  ■ 


and 


(9a) 


(9b) 


Substituting  (8)  into  (6)  we  get 


(10a) 


~~  -ji>  + 


cK/J 


e 


fltxi 


Upon  a  little  rearrangement  these  equations  take  on  an  identical  form: 


and 
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Because  of  the  properties  of  the  operators 


this  equation  can  be  written  in  a  more  compact  form: 


in,n 


K 


(12) 


where  m,  n  now  run  over  all  states  and  over  all  particle  numbers,,, 
and 


(13) 


The  fact  that,  in  the  Ep-rep,  the  "retarded"  and  "advanced" 
temperature  Green's  functions  are  identical  is  a  reminder  that  the 
analogy  between  the  temperature  variable  and  time  is  not  exact. 
This  will  become  clear  when  we  discuss  the  retarded  and  advanced 
time  Green' s  functions  in  the  next  section. 
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(c)  The  Ep-rep  of  the  Time  Green’s  Function  for  Finite  Temperatures 
We  now  summarize  the  salient  features  of  the  p-rep  for 
the  time  Green's  functions.  Proceeding  in  much  the  same  way  as 
was  done  for  the  temperature  variable  in  the  previous  section  we 

get 

itH  -itH 

f(ne 


itH  y  _//// 
y  (x)  -  e  f  c r )  e 


(14a) 


(14b) 


The  Green's  function  in  the  momentum  representation  is 


U’t-f) 


K  ^  ,  fJ aft 


(15a) 


6  Djadi  forA<0.  (15b) 


ti 


a  cl 


Comparing  equations  (12)  with  (15)  we  see  immediately  that  G  (  /  ) 

K 

(  y  ^  Q  )  is  the  analytic  continuation  of  Or  (if)  (if2o):0+-tf 
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(cf.  Fradkin,  1959).  Thus,  given  GK(  V  ),  we  can  immediately 
write  down  the  corresponding  G  (  ft  ).  This  is  not  surprising  in 
view  of  the  parallel  formulations  which  we  have  presented  in 
Chapter  III  for  the  temperature  and  time  Green' s  functions. 


We  now  define  the  following  functions: 


(a) 

G  K  U) 


‘  G.  I/)  *>°  > 


0 


/ 


i  4  0  • 


(16) 


Cr  u  ft) 


0 


ft  >  o  > 


;  GKfr)  /(O'. 


a. 


,tr) 


The  function  G  K  (ft  ),  respectively  G  K  (ft-  ),  can  be  readily 
expanded  into  its  Fourier  decomposition  in  the  lower*  respectively 
upper*  half  of  the  ft  plane;  the  result  is 


/ 


G  (?)  =  v7  e'  —  „ 


ft,  >%,Z<0,W »> 


=  ^4  e 

//a  a 


"4  ^4  f 


c  _  F  '  -  7  //aa 

CH,  *  *  A 


3,.  y3ml>0,{\ 7b) 


where 


-  oo 


&la'r\z)  =  |  Gft'Wexftiftd/  , 


OO 
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and  where  the  notation  of  the  previous  section  has  been  employed. 

The  singularities  of  these  expressions,  when  regarded  as 
functions  of  the  complex  variable  j?  ,  occur  only  along  the  real 
axis.  Furthermore  in  the  limit  of  an  infinitely  large  system,  as 
is  considered  here,  each  point  is  usually  a  limit  of  poles  and, 

therefore,  the  real  axis  is  in  general  a  branch  line  (cf.  Luttinger,  1961). 

{Q.)  (r  \ 

Hence  G  ' f.  (  t  ),  respectively  G  K  (  %  )  ,  can  be 

analytically  continued  into  the  whole  lower,  respectively  upper,  plane. 


(59) 


(d)  Comparison  Between  the  Ep-rep  of  the  Temperature  and  Time 
Green's  Functions  at  Finite  Temperatures. 


function 

&/<?,) 


We  now  split  the  Ep-rep  of  the  temperature  Green 
(equation  (11)  )  into  two  parts :  GK  (  )  = 

+  itA)  ;  where 

f\ 


8 


(18a) 


and 


fc  (ff)  -  D, 

^  a/  1 

A  a  ft 


(18b) 


In  the  above  expressions  we  have  used  equation  (9a)  and  the  assumption 

that  the  summations  over  a  and  a  run  over  the  same  quantum 

states.  The  last  assumption  is  based  on  the  fact  that,  physically,  we 
expect  the  /\l  particle  and  fj  t  /  particle  systems  to  be  characterized 
by  the  same  quantum  numbers.  As  in  the  case  of  the  Ep-rep  of  the 
time  Green's  functions,  these  expressions  (regarded  as  functions  of  a 

complex  variable  T  <  te  —  *  >  >  are  analytU  ever>rwhere 

except  on  the  real  axis. 

Comparing  (18a)  with  (17a)  and  likewise  (18b)  with  (17b), 
observe  that  these  expressions  are  simply  related  in  the  upper. 


we 
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respectively  lower,  planes.  From  these  observations  we  can  readily 
see  that  i  •  t)  represents  the  Ep-rep  of  the  advanced 

Green's  function  in  the  upper  plane,  and  -  <  &h  it)  the 

retarded  Green's  function  in  the  lower  plane.  (The  factor  l 
arises  from  the  definition  (16)  ).  It  is  difficult,  and  perhaps 
impossible,  to  extend  the  identification  to  the  whole  plane  at  finite 
temperatures  because,  in  this  case,  the  whole  real  axis  could  be  a 
branch  line  in  the  complex  plane.  The  latter  is  due  to  the  presence, 
in  general,  of  continuous  spectra  of  energies  and  of  degeneracies  at 
the  various  excited  states.  However,  in  the  zero  temperature  limit 
one  can  extend  the  identification  to  the  whole  complex  plane.  This 
is  discussed  in  the  next  section. 
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(e)  Ep-rep  of  the  Green's  Functions  in  the  Limit  of  Zero  Temperature. 
In  the  limit  of  zero  temperature,  equations  (15)  reduce  to 


e 


■i/far  bni 


3, 


K 


o,n 


where  we  consider  an  v  particle  system  whose  ground  state  is 
given  by  (  W,  and  its  ground  state  energy  by  /TW)0  .  The 

sum  over  n  runs  over  all  states  and  all  particle  numbers;  n  =  0 
denotes  the  ground  state  which  is  assumed  to  be  non -degenerate . 
That  equations  (19)  are  correct  is  clear  inasmuch  as  all  other  terms 
in  (15)  are  exponentially  smaller  in  the  p  — ^  OO  limit  as 

compared  to  the  terms  given  by  (19).  Proceeding  as  in  section  (c) 
we  get 

=  (M flK -JTf—J  %  I ^  ■  f‘r  W  ?  ° 


and 
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for  °)<n  £  <  0 


(20b) 


The  zero  temperature  limits  for  the  Mtemperature"Green '  s 
functions  in  the  Ep=rep  (cf.  Luttinger,  1961)are  given  by 


/ 

h 1 


(21a) 


and 


/ 


L,.  -tf'T 


(21b) 


With  regard  to  the  derivation  of  equations  (21),  we  merely  remark 
here  that  in  the  limit  j2  — Oo  ,  ^  +/** 

becomes  a  continuous  variable,  and  hence  equation  (21b)  becomes  an 
integral  transform  equation.  The  details  of  the  derivation  are  given 
in  Luttinger  (1961)  and  in  Luttinger  and  Ward  (I960). 
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If  we  regard  Q  '  (  $  )  as  a  function  of  the  complex 

variable  ^  we  note  i™ediately  that  the  fanction  is  analytiC 
everywhere  except  on  the  real  axis.  But  in  the  case  of  interest 

nQw  (i.e.  $  OO  )  we  shall  show  that  both  functions  are 

analytic  also  on  infinite  segments  of  the  real  axis.  Physically 

this  can  be  seen  as  follows:  -  the  chemical  potential  for  an  f\j 

particle  system  can  be  interpreted  as  the  change  in  the  energy  of 
the  system  upon  removal  of  one  particle  at  fixed  3  and 

volume:  /< N  =  Ac  '  F«-i>o  '  where  both  energies  are 

ground  state  energies.  Furthermore,  for  large  A  we  expect 

(e.g.  Pines,  1961)  that  A  */*»+'  ~  A*''  »  which  we 

denote  simply  by  /<  •  The  correction  to  this  equation  is  of 

the  order  of  '//V  and,  in  our  case,  can  be  neglected. 

With  this  equality  we  can  write  the  energy  differences  occurring 
in  equations  (21)  and  (20)  as  follows  (cf.  Pines,  1961): 


^ H,  D 


E .  -  E„t% 


E  -  -b  (Ejh.q  1  ?  > 
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and  thua  equations  (21a)  and  (21b)  are  regular  on  the  real  axis 
in  the  region  —  oo  to  /<  (  jM  to  -t  oo  ).  This  clearly 

holds  for  equations  (20)  as  well.  Hence  we  conclude  that  in  the 
zero  temperature  limit,  iG  *  (  ^  )  represents  the  Ep-rep  of  the 
advanced  Green's  function  in  the  whole  complex  ^  plane,  and 

similarly  -  i  Q-  (  ^  )  represents  the  retarded  Green's 

K 

function.  The  reason  that  this  identification  is  possible  in  the 
zero  temperature  limit  is,  clearly,  the  possibility  of  analytic 


continuation  of  G 


(a.,r> 

K 


(z)  into  the  upper  and  lower  planes,  since 


in  this  limit  (  p  oo  )  the  real  axis  is  not  a  branch  line  from 


-  oo  to  -f  OO  . 
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(f)  Spectral  Representation  for  the  Green's  Function  in  the  Zero 
Temperature  Limit. 


For  an  infinite  system  the  energy  spectrum  is  expected 
to  form  a  continuum,  i.e.  the  poles  of  equation  (25)  come  close 
together  and  the  real  axis  is,  in  general,  a  branch  line  (luttinger, 
1961).  For  the  case  of  very  closely  spaced  energy  levels  it  is 
appropriate  to  convert  the  summations  into  integrals.  This  is  done 
customarily  (e.g.  Pines,  (1961),  Luttinger,  (1961))  by  introducing 
the  spectral  functions 


Mk,*)) 

BM 


Here  and  c  are  defined  by 

equations  (9)  and  (13).  With  these  definitions  the  formulae  (25) 
can  be  written  as 


(f) 


r 


A 


/\  (k,co)  d^j 

60  -  f 


-  oc 


) 


H(k^)  d co 


(a)  -  f 
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where  60  plays  the  role  of  the  energy  differences  in 
equations  (25).  The  association  of  the  time  and  temperature 
Ep-reps  will,  clearly,  go  through  in  much  the  same  way  as 
before . 
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CHAPTER  V:  CONCLUSIONS  AND  OUTLOOK 

We  begin  this  chapter  with  a  short  historical 
survey  of  the  subject.  This  survey,  which  one  might  normally 
have  expected  in  the  introductory  chapter,  has  been  deferred  for 
presentation  here  because  the  literature  in  this  field  is  relatively 
recent  and  the  problem  of  discriminating  between  the  contributions 
of  various  authors  is  more  simply  handled  after  the  general  methods 
have  been  discussed. 


The  temperature -dependent  many  body  theory 
aims,  as  does  its  zero  temperature  counterpart,  to  obtain 
expressions  for  the  energy  spectrum  and  the  mean  density  of  a 
system  when  the  mutual  interaction  of  the  particles  is  given. 

The  T  =  0  case  was  solved,  in  a  certain 
sense,  by  Brueckner  and  his  collaborators  (references  are  given  in 
Goldstone's  paper);  however,  their  approach  was  ambiguous  and 
difficult.  Goldstone  (1957)  was  the  first  to  attack  the  problem  using 

the  time -dependent  techniques  associated  with  Feynman  diagrams 

> 

borrowed  from  field  theory.  His  approach  was  further  simplified 
by  Bloch  (1958).  The  modern  approach,  wherein  extensive  use  is 
made  of  the  Green' s  function  formalism,  was  initiated  by  Klein  and 
Prange  (1958);  these  authors  extended  their  considerations  to  the 
T  0  case. 
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The  theory  of  finite  temperature  MBP  is  intrinsically  more  complex 
than  the  T  =  0  case  because  of  the  inherent  requirement  for 
statistical  averaging;  the  theory  cannot,  at  present,  point  to  a 
success  comparable  to  the  theory  of  Brueckner  and  Goldstone  for 
zero  temperature .  Nonetheless,  some  important  achievements 
have  been  made:  (a)  It  is  now  possible  to  obtain,  in  the  T  — >-  0 
limit,  the  Goldstone  expression  for  the  ground  state,  without 
recourse  to  time -dependent  theory  (the  Green's  function  formalism 
can  be  considered  here  as  a  concise  formulation  for  the  time- 
dependent  approach).  This  point  was  clarified  by  Luttinger  and 
Ward  (I960),  who  considered  a  perturbation  expansion  wherein 
a  temperature -like  variable  (denoted  by  u  in  this  thesis) played 
the  role  of  time,  (b)  The  semi -phenomenological  theory  for  Fermi 
liquids  proposed  by  Landau  (1957)  was  put  on  a  firm  theoretical 
footing  by  Klein  (1961b),  who  used  the  "time " -Green '  s  function. 

This  theory  was  also  derived  from  first  principles  by  Balian, 

Bloch  and  DeDominicis  (1961)  who  used  the  "generalized  Wick 
contractions"  first  devised  by  Bloch  and  DeDominicis  in  19  58; 
this  method  also  employs  the  U  variable  explicitly,  (c)  Luttinger 
(I960)  obtained  a  temperature -dependent  extension  for  the  Hartree-Fock 
approximation  which  reduces  (see  qualifying  remarks  in  Tnouless, 

1961)  to  the  standard  Hartree-Fock  expression  in  the  zero  temperature 
limit.  This  approach  once  again  uses  the  u  variable  as  an 
auxiliary  tool  for  the  discussion. 
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AGD  (1959)  and  Fradkin  (1959)  re-examined  an  approach  to 
the  T-dependent  MBP  which  had  been  suggested  originally  by 
Matzubara  (1955).  These  authors  extended  Matzubara's 
u -variable  approach  and  pointed  out  an  interesting  parallelism 
between  the  rdles  played  by  the  time  and  the  u-variable  in  the 
theory.  For  example,  AGD  considered  the  retarded  "time 
Green1  s  function  and  showed  it  to  be  the  analytic  continuation 
of  the  "temperature"  Green's  function.  Fradkin  in  a  relativistic 
formulation  of  the  problem  considered  the  "temperature" Green '  s 
function  in  some  detail  and  outlined  a  scheme  for  the  determination 
of  the  energy  spectrum  of  a  system  in  terms  of  it. 


As  may  be  clear  to  the  reader,  the  main  interest  in  this 
thesis,  apart  from  an  elaboration  of  the  Green's  function  itself, 
has  been  the  study  of  the  curious  and  suggestive  parallelism  between 
the  roles  played  by  the  time  and  the  u-variable  in  the  T-dependent 
MBP.  We  outlined  the  conventional  theory  of  the  Green's  function 
(i.e.  the  "time"Gr een '  s  function  )  as  it  has  recently  been  formulated 
for  the  MBP.  Then  a  new  variable,  T  ,  was  introduced  into 
the  theory;  T  was  shown  to  play  the  two  roles,  viz.  time  and  u 
in  the  conventional  theory.  The  purpose  of  this  formalism  was  to 
make  the  parallelism  between  the  u  and  t  variables  manifest, 
but  it  also  enabled  us  to  rederive  results  which  had  been  hitherto 
obtained  by  Luttinger  and  Ward  (without  recourse  to  the  Green  s 
function  formalism).  It  is  at  once  evident  that  our  derivation  is 
simpler  and  also  appreciably  shorter. 
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The  striking  parallelism  between  the  roles  of  the  time  and 
the  temperature  variable  in  the  T~  formalism  (and,  of  course, 
in  general)  is,  as  one  would  expect,  not  complete:  for  one  thing 
the  u -variable  occurswith  an  imaginary  argument  in  T  while 
the  time  occurs  with  a  real  argument;  also  there  is  the  fact  that 
the  variable  u  is  confined  to  a  finite  region  (  1 LL l  <  P/fc  ) 

while  t  is  unbounded.  These  differences  are  manifested  most 
clearly  in  the  Ep-rep  where  we  showed  that  while  one  is  able  in 
the  usual  way  to  distinguish  between  the  "retarded”  and  "advanced' 
time  Green's  function*,  this  distinction  is  not  present  for  the 
"temperaturel,Green' s  functions.  However,  we  have  shown  that 
in  the  Ep-rep  the  "temperature"Green '  s  function  can  be  split  into 
two  parts,  one  of  which  can  be  associated  most  naturally  with  the 
advanced  and  the  other  with  the  retarded  time  Green's  function  . 

It  is  our  hope  that  the  ~T  formalism  can  be  used  to 
investigate  a  fur  ther  -going  connection  between  the  t  and  u 
variables.  We  have  also  speculated  briefly  on  the  possibility  that 
the  introduction  of  non -diagonal  elements  in  the  resultant  Green' s 
function  might  correspond  to  terms  describing  phenomena 
associated  with  the  time  development  of  non-equilibrium  conditions*’, 


*  e.g.  Professor  Y.  Nambu  in  a  private  discussion  suggested  that 
insertion  of  some  appropriately  chosen,  non -diagonal  terms  info 
our  Green's  function  might  lead  to  a  new  formulation  of  the  problem 


of  thermal  conductivity. 
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however,  our  predilection  lies  in  farther  investigation  of  the 
above-mentioned  analogy  between  time  and  the  temperature 
variable  u.  In  particular,  in  continuing  the  investigation,  the 
splitting  of  the  temperature  Green's  function  in  the  Ep-rep  into 
"retarded"  and  "advanced"  parts  will  be  examined  in  other 
representations  to  see  whether  a  deeper  significance  to  the 
analogy  between  the  u  variable  and  time  can  be  uncovered. 
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APPENDIX  A:  VERIFICATIONS  OF  SOME  FORMULAE  GIVEN 

IN  CHAPTER  I  I 

The  proofs  and  discussions  given  below  lean  heavily 
on  the  lucid  review  article  by  Kato  et.  al.  ,  (I960). 

(i)  Source  terms  in  the  Lagrangian  density. 

Let  us  consider  a  closed  system  with  no  external 
sources.  Since  the  Lagrangian  will  then  have  no  explicit  dependence 
on  the  time,  the  resultant  Hamiltonian  will  possess  a  spectrum  of 
stationary  eigenvectors  which  we  label  by  ^  .  We  now  add 

the  following  terms  to  this  Lagrangian, 

-  J  +U)  ^  (*)  -  ]  Cx)  ) 


Source 

where  J(x)  is  an  infinitesimal  external^ which  is  independent  of  time 
for  t  «  t  and  for  t  t  .  As  an  examples  consider 

J(x)  to  have  a  step  function  time-dependence:  zero,  say,  for 
-  oo  <  t  <.  t*  ,  non-zero  but  infinitesimal  for 

t  <C  t  <C  00  •  The  Hamiltonian  resulting  from  this 

modified  Lagrangian  will  not  have  a  stationary  spectrum  of 

,  i 

eigenvectors.  Nonetheless,  we  expect  that  for  t  <3C  *  it 


does  have  the  spectrum 


j  ITU 


.  I 
I 


while  for 


t  »  t 


it  has 


a  new  spectrum  which  we  label  |  »  further,  in  view  of  the 

source  being  infinitesimal,  we  expect  and  assume,  that  there  exists 
a  one  to  one  correspondence  between  the  spectra  ^  til-  ^ 


and 


A  2 


|  IT) j.  j  such  that,  as  we  set  the  source  equal  to  zero,  there 

exists  an  which  becomes  identical  with  IT?  •  for  every  tn- 

in  the  spectrum.  More  succinctly,  we  demand  for  physical 
Hamiltonians : 


5  V 


IT?; 


3  m  m-  as  J 


o 


(Al) 


(ii)  Summation  over  initial  and  final  spectra  (equation  (8)  ): 

With  the  correspondence  (Al)  of  to  fT?  •  for 

every  In-  ,  the  meaning  of  the  summation  defining  equation  (8b) 
can  be  made  explicit  by  writing 


P  H 


=  j  m-^>  60 


A. 

m; 


m.  \mf 


(A  Z) 


where  {jj  is  the  weight  factor  specified  by  equation  (8). 

1 1 

(iii)  The  functional  derivative  :  &  Th  /  b  J(x)  . 

Let  us  now  vary  the  Lagranqian  by  letting 
'Jlx)  — J  (x)  -+  £  >  where  it  is  convenient  to  assume  that  £ 

is  such  as  to  make  £f(x)  an  infinitesimal  of  higher  order  than  J(x), 

and  f(x)  is  an  arbitrary  function  in  the  sense  of  equation  (1);  for 

.1 

convenience  we  assume  it  to  be  localized  in  time  around  t  .  We 

label  the  independent  variations  of  J(x)  and  j\x)  by  S  and  b>  t  , 
respectively;  e.g.  X  reads:  the  variation  in  X  due  to 
the  variation  in  Jf(x).  With  these  preliminary  remarks  we  can  take 


A3 


the  functional  derivative  of  Tr.  jO  with  respect  to  the  source 
functions  as  follows: 

b Tr.  p  =  odm.  <mf  |  m;>  , 

hi; 


which,  through  Schwinger's  principles  becomes 


5 


r. 


F 


60, 


—  -  t  /  ^  m 

m  ■ 


)  d  A'  d  6a)  £  }(*)  j  177 .  (A  3) 


On  the  other  handi 


(5>T r.p-  i  r 0 


Tr.  0 


j*£f  I 


J 


f  k>  Tr.  f>  ^  . 

— —  £f(*)d* 
d  J  W 


( A4) 


in  the  limit  £  0.  The  last  identity  follows  from  the  definition  of 

functional  differentiation  (equation  (1)  ). 

Comparing  (A3)  with  (A4)  we  get,  since  f(x)  is  arbitrary, 


=  -i  J  com,  FmiW(x)\ m-)  • 

6T<x)  Ft 


(A  5) 


A4 


(iv)  The  functional  derivative: 

Since  the  source  functions  must  be  spinors  for  the 
Lagrangian  to  be  scalar*  they  anticommute  with  other  spinors  such 
as  Q  (x).  Hence,  in  considering  the  matrix  element  j  Q(*)  j 

we  get  two  distinct  possibilities  depending  on  whether  the  time 
associated  with  Q  (labelled  by  t  )  is  earlier  or  later  than  the  time 
of  variation  of  the  source  (t1  ).  For  t  )>  t '  we  choose  a  time  t, 
such  that  t  >  t(  >  t'  ,  so  that  at  t,  we  have  once  again  a 

spectrum  satisfying  the  stipulation  (Al).  We  label  this  spectrum 
by  \  m(t()  ^  I  hence*  we  have  immediately: 

6  j  (*) 

where  plays  a  role  of  fo<>  in  Schwinger's  dynamical 

principle.  The  same  discussion  leading  to  (A5)  now  leads  to: 


£ 

b  ft  ft 


•(mfl  Q(x)im,}  ~-i(milQ(x)^ii;x)\miy  ■ 


(A6a) 


A5 


An  analogous  procedure  for 


t'  >  t 


leads  to: 


& 

STU) 


(nyl  Q(x)  |  m^> 


m 


fU)  QC*) 


m. 


(A  6b) 


where  the  sign  change  is  due  to  the  anticommutativity  of  the  source 
with  the  spinor  Q(x).  Equations  (A6)  can  be  written  as: 


-i(nylT  QWfU* 


IV; 


(A  7) 


where  T  is  the  chronological  operator  of  Wick. 


(v)  The  statistical  Green's  function. 

The  Green's  function  of  our  problem  was  defined  by: 


(y-(x,x)  — 


8<f> 

8J(x) 


J  ~o 


where  the  expectation  values,  signified  by  the  angular  bracket,  is 
defined  in  equation  (8): 


( 


Hence> 


b 


& 


kV/  SJM)  lr/° 


r. 


r 


which  upon  application  of  (A7)  and  (A3)  reads: 


blU) 


where  now 


f 


(x)  plays  the  role  of  Q(x). 


If  we  now  take  the  limit 


tTl fbTl i  (by  letting  J  0),  the  second  term  vanishes  and 

we  get  finally: 


in  agreement  with  equation  (9). 


A7 


(vi)  The  Schwinger  identity. 

The  Schwinger  identity  states  the  following  operator 


relation: 


t  * 


/ 


The  proof  can  be  obtained  in  the  following  way: 


where  £)  (x)  is  the  usual  step  function 

I  x  >  o 

9(x)  = 

o  X  <  0 


Differentiating  (A8)  partially  with  respect  to  the  time  t  t  we 
get: 


t  ^  \x)f[x)  ~  Qli~t)  y  (X)  j(x)  y 


(AS) 


(A9) 


A8 


where  we  have  used  the  well-known  formula:  o> x  Bi-*)  -  o  (.  *  )  . 

Using  the  commutation  relation  of  ~\J/  and  y  f  we  can 

combine  the  first  and  third  terms  in  ( A 9) ,  while  the  second  and 
fourth  terms  can  be  combined  using  the  definition  of  the  chronological 
operator  T.  These  combinations  result  in  the  Schwinger  identity 
(equation  (13)  ) . 
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APPENDIX  B.  COMPLETENESS  OF  THE  SET 


d> 


i 


The  set  0  was  defined  by  (equation  (20), 

'  >  a 


Chapter  III): 


(*>o)  -  ~7  -<*/>  <?-r 


The  spatial  factor  is  the  familiar  Fourier  one,  and  as  such,  requires 
no  further  comment.  The  factor  containing  ij  and  involving 

TC L  + possesses  certain  symmetry  properties; 


->  f  ~ 


V  P 

hence,  it  is  necessary  to  investigate  whether  this  set  is  sufficiently 
complete  for  our  purposes.  We  consider  the  spatial  Fourier 
component  of  (j-  (x>X)  and  show  that  it  satisfies  the 

requirements  which  warrant  an  expansion  in  terms  of 


8 


1  **p[-  5  *]  - 


We  had, 


&j>(u-,U)  = 


T 


a 


/ 


V 


u  -  u. 


): 


Written  more  fully  this  equation  reads  ( 
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_  -ph'i-  o(Ai  _  2///  f 

(r  M  =  -  J)  !r.  6  6  ^  £  dp  {or  J<  o  ,  (bu) 


/-  /  ,  -r.  -7-  -0+o(/J  -jH  f  -M 

=  Jj  !n  e  6  dp  C  Cip  {or  7) >0  ,  (Bib) 

(see  equations  (3)s  (4)  of  Chapter  IV). 

Now  consider  the  expression  (Luttinger*  1961;  Bloch*  1960); 


This  expression  has  the  following  properties: 

SpM  -  %  (V-y3J  =^(1*)  jor  O<0,i)>0  , 

Sfih  -  ~  f°r  V'>0,  v<0  . 


All 


These  properties  can  be  established  from  equations  (Bl),  by 
using  the  cyclic  invariance  of  the  trace  and  the  following  properties 
oj:  the  operators  Cl ^  and 


aK(  a/-i)  =mqk 
a+K  ( h + 1 )  =  dciK 


This  symmetry  property  of  Q  (  ~5  )  is  summarized  by  the 

P 

statement:  Cj  (  id  )  is  periodic  with  a  period  in  the 

interval  (  —  }  ^  )■  Since  only  this  interval|enter s  into  our 

calculation,  we  can  assume  0  (  V  )  to  be  periodic  everywhere 

P 

This  observation  implies  that  G  (  1)  )  can  be  expanded  in  the 

Jf> 

Fourier  series 


and  hence 


Clearly  a  similar  proof  will  go  through  for  Q-  (  V  )»  the 

non-interacting  Green's  function. 


Now  equation  (14b)  of  Chapter  III  can  be  expressed  by: 


which,  in  turn,  can  be  written  in  terms  of  the  inverse  Green's 
function: 


Naturally,  if 


can  be  expressed  in  terms 


